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' Abstract: We study the constraints on the inflationary parameter space derived from 

^ ■ the 3 year WMAP dataset using "slow roll reconstruction" , using the SDSS galaxy power 

' spectrum to gain further leverage where appropriate. This approach inserts the inflationary 

Q^. slow roll parameters directly into a Monte Carlo Markov chain estimate of the cosmological 

Q , parameters, and uses the inflationary flow hierarchy to compute the parameters' scale- 

|H \ dependence. We work with the first three parameters (e, r] and ^) and pay close attention 

^ \ to the possibility that the 3 year WMAP dataset contains evidence for a "running" spectral 
index, which is dominated by the ^ term. Mirroring the WMAP team's analysis we find 
that the permitted distribution of ^ is broad, and centered away from zero. However, when 

j_j ■ we require that inflationary parameters yield at least 30 additional e-folds of inflation after 



the largest observable scales leave the horizon, the bounds on ^ tighten dramatically. We 
make use of the absence of an explicit pivot scale in the slow roll reconstruction formalism 
to determine the dependence of the computed parameter distributions on the pivot. We 
show that the choice of pivot has a significant effect on the inferred constraints on the 
inflationary variables, and the spectral index and running derived from them. Finally, we 
argue that the next round of cosmological data can be expected to place very stringent 
constraints on the region of parameter space open to single field models of slow roll inflation. 
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1. Introduction 

The primordial perturbation spectrum poses two distinct problems for cosmology. Obser- 
vationally, we wish to extract the spectrum's amplitude and scale dependence from data, 
while theoretically we seek to understand the origin of the perturbations. After the release 
of the 3-year WMAP dataset [||, ^ ^ there is considerable optimism that we will soon mea- 
sure the departure from scale invariance in the spectrum at significance levels greater than 
3a. ^ On the theoretical front, inflation ^, |6| remains the leading theoretical paradigm 
for understanding the very early universe and the origin of the perturbations. In general 
terms, the latest round of cosmological data has confirmed the broad predictions of infla- 
tion - the universe is very close to being flat, it has a roughly scale-invariant perturbation 
spectrum and {TE), the cross correlation between the E mode of the polarization and the 
temperature, has the form predicted by inflation. 

Beyond this rosy picture, however, lies a number of qualifications. Firstly, the observed 
form of {TE) is not a unique prediction of inflation, but is also a generic feature of models 
where the perturbations are laid down during a collapsing phase prior to a bounce, such 
as the ekpyrotic scenario 0].^ Secondly, inflation can be implemented in a vast number of 
different ways, and each mechanism produces a potentially distinctive perturbation spec- 
trum. Consequently, the inflationary paradigm does not make an unambiguous prediction 
for the detailed form of the perturbation spectrum. Some authors have posited that in- 
flation has a "default" perturbation spectrum [Q]. The underlying arguments leading to 
this position rely on simplicity and naturalness, rather than a watertight "no go" theorem. 

^WMAP is an acronym for the Wilkinson Microwave Anisotropy Probe. Here WMAPI refers to the first 
year's data from the satelhte, and WMAPII refers to the second WMAP data release, which is drawn from 
the first three years' data. 

^The (TE) signal is the same for any mechanism that induces adiabatic perturbations that are correlated 
on scales much larger than the post-recombination horizon volume. "Bouncing" models do this by laying 
down perturbations during the contracting phase, when the comoving horizon size is shrinking. 
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Consequently, inflation itself will not be ruled out if the perturbation spectrum turns out 
to differ from this simple expectation, but it is certainly true that many well-known models 
would be excluded. 

Looking at data, it is always tempting to focus on possible anomalies - parameters 
that differ noticeably from their "expected" values, but for which the discrepancy falls well 
short of a 3a detection. Since the release of the original WMAPI dataset 0, three 
aspects of the power spectrum have received considerable attention: the "low" values of 
Ci at small i, the glitches around the first Doppler peak and at £ ~ 40, and the apparent 
evidence for a running, or scale-dependent, spectral index that arose when the WMAPI 
data was combined with other probes of the power spectrum at shorter scales. The status 
of all of these anomalies has evolved since the WMAPI release. The quadrupole is still 
~ 2a lower than the LCDM best-fit value, but the octupole has moved closer to the LCDM 
value, thanks to an optimal extraction of the spectrum. The glitch near the first Doppler 
peak is not present in the WMAPII dataset with its smaller errorbars, while the ^ ~ 40 
discrepancy is still seen |12, Finally, the WMAPII dataset contains a preference for a 
running spectral index when viewed alone, in combination with other CMB data, and/or 
large scale structure information such as 2dFGRS |14] and the Sloan Digital Sky Survey 
|15]. However, a reanalysis by Seljak et al., which includes SNla data and information 
from the Lyman-a forest, puts a tight constraint on any running ||T^ — this result seems 
to hinge on the inclusion of the Lyman-a data. Further analysis of a possible running index 
is provided by |17]. 

In this paper, we focus on the running of the spectral index, and the constraints it 
places on the inflationary parameter space. Conventionally, one specifies the spectrum at 
a pivot point in terms of the spectral index and the running Inflation predicts Ug 

and in terms of the slow roll parameters, which are functions of the potential and its 



derivatives. In ||1^ we showed that if the running is as large as the central value extracted 
from the WMAPII dataset, no simple model of inflation can simultaneously provide both 
this running and a sufficient amount of inflation.^ Currently, there are simple models of 
inflation within the errorbars, but if the constraints from the data continue to tighten 
around the median value of the running, they will be ruled out. We arrived at this result 
by considering the Hubble Slow Roll hierarchy, otherwise known as the inflationary flow 
equations 2C, 21, 22|. These codify the scale dependence of the slow roll parameters: 
given the values of the slow roll parameters at a single moment during the inflationary era, 
we can compute their values at all other times. A further virtue of this approach is that 
when the infinite slow roll hierarchy is truncated at finite order, the truncation is preserved 
by the evolution equations. Using the slow roll hierarchy means that we effectively dispense 
with the pivot point, other than when comparing our results to treatments using the 
standard Ug and parameterization of the spectrum. 

As we showed in |22], we can include the slow roll parameters iii the cos- 

mological parameter estimation process, and entirely avoid the use of the spectral index 
and its running. Further, the slow roll hierarchy can be solved to produce the underly- 

^Which we take to be a minimum of 30 e-foldings, and by "simple" we mean a single inflaton, minimally 
coupled to gravity, whose potential can be well specified in terms of its first three derivatives. 
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ing inflationary potential, leading us to dub this approach "slow roll reconstruction."^ Of 
course, this presupposes that the perturbations were produced by inflation, so slow roll 
reconstruction does not replace fits to n<j and -j^^, which are purely empirical charac- 
terizations of the spectrum and contains no strong theoretical priors.^ However, if we do 
adopt an inflationary prior, fitting to the slow roll parameters directly and using the flow 
equations to obtain their scale dependence ensures that we gain the maximum possible 
leverage from imposing this constraint. It is well known that the scalar and tensor pertur- 
bations obey a consistency condition at lowest order in slow roll,^ but there are further, 
more subtle correlations at second order and beyond which tie together the running and 
the scalar spectral index [^. By using slow roll reconstruction we automatically include 
these correlations in our analysis. Moreover, this removes a potential source of ambiguity 
inherent in the introduction of a scale dependent spectral index - which is how to best 
parametrize the scale dependence? While is arguably the simplest measure of the 
scale dependence, it is by no means unique. In particular, by adopting it we are implicitly 
assuming that while varies with k, -^j^ does not. By specifying the spectrum in terms 
of the slow roll parameters we ensure our spectrum has the scale dependence predicted by 
inflation, and the primary limitation is the highest order slow roll parameter we include 
in our fit. While this is obviously a choice that will be guided by the quality of the data, 
it means that we can describe the scale dependence in a relatively unambiguous way.^ Fi- 
nally, since we obtain the overall form of the potential as a by-product of this analysis, we 
can include constraints based on the total number of e-folds of inflation. 

The closest analogue to slow roll reconstruction is the approach of Leach and collabo- 



rators using the Horizon Flow Functions 27, who compute the spectral index and 



running as functions of the slow roll parameters. As in slow roll reconstruction, one has an 
explicit inflationary prior, and the slow roll parameters appear directly in the parameter 
set. The advantage of this approach is that it can be implemented more simply in existing 
codes, since it translates between the slow roll formalism and and -j^^, whereas in 
order to implement slow roll reconstruction we have to numerically solve the differential 
equation that maps the field evolution into the comoving wavenumber k.^ However, slow 
roll reconstruction can be more easily extended to include higher order terms in slow roll, 
and copes naturally with heterogeneous datasets that span a wide range of wavelengths - 
this is not impossible in the formalism of p6| , p^ , 28 1, but would require more algebraic 

*The overall amplitude of the spectrum is fixed by an arbitrary constant of integration, which can only 
be measured directly if the tensor spectrum is observed. 

^The index and running are effectively the first two terms in a Taylor series expansion, so their use 
does make the tacit assumptions that the spectrum is a smooth, differentiable function of the comoving 
wavenumber k. 

^In practice, almost all estimates of the cosmological parameters assume that the amplitude and slope 
of the tensor spectrum have the correlation predicted by minimal models of infiation. Since the tensor 
spectrum is merely bounded from above, rather than observed directly, this constraint makes a Monte 
Carlo Markov Chain analysis more tractable without excluding viable models. 

'^This issue mirrors the concern when the dark energy equation of state is allowed to vary with time, as 
there is no unique way of specifying that time dependence | p^ , p5| . 

®In practice this is not computationally expensive, especially when compared to the numerical evaluation 
of the CMB spectrum and likelihood codes. 
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labor. Secondly, slow roll reconstruction can naturally include constraints based on the 
total duration of inflation. Finally, and perhaps most importantly, because 27, 28| 
use n<j and they are effectively Taylor expanding the power spectrum round a fixed 

pivot scale, /cq- Conversely, we can easily post-process our chains to present parameter 
constraints at an arbitrary pivot. In p9| ], the formalism of |26, 28 1 is applied to the 



WMAPII dataset to produce constraints on the slow roll parameters, and this work is 
directly comparable to the results we present below. Slow roll reconstruction also draws 
on the stochastic Monte Carlo reconstruction algorithm, described in ||2^ . This procedure 
does not yield unambiguous estimates of the slow roll variables, but it permits estimates 



of the impact of higher order terms in the slow roll expansion |22| and has been applied 



to cosmological data in jl^, |3y, |3^, a similar approach can also be found in 



An alternative approach of computing exact power spectra for specific potentials has been 



apphed to the WMAPII data by |3§. 



Our specific goal here is to extend our analysis in ||2^ to include the third slow- 
roll parameter, ^, in a Monte Carlo Markov Chain [MCMC] analysis. This parameter 
dominates the running of the spectrum, so it must be incorporated in any discussion of 
a running spectral index, especially in the light of the hints in the WMAPII data that 
the spectrum has a non-trivial scale-dependence. We will also include constraints based 
on A^, the number of e- folds of inflation remaining after the fiducial scale has left the 
horizon. We consider three different sets of chains: WMAPII + > 30, WMAPII + 
SDSS + A > 30, and WMAPII + SDSS with no N constraint. While we put a cut on 
N, we do not take the further step of correlating A^ with the inflationary energy scale. A 
low N generally implies that inflation occurs well below the GUT scale but making this 
statement precise requires introducing assumptions about the physics of reheating, which 
we wish to avoid. ^ Roughly speaking, however, imposing A > 30 is equivalent to setting 
the inflationary scale at or above the TeV scale. This ensures that inflationary energy 
scales are not probed by precision electroweak measurements, and allows baryogenesis to 
be driven by electroweak physics. Given that the GUT scale is often regarded as the 
"natural" inflationary scale, assuming that inflation happens above the TeV scale is a very 
mild theoretical prior. However, demanding A > 30 leads to very tight bounds on ^. 
Translating our results into bounds on yields a tighter constraint than that of Seljak 



et al. 1 16] when we add the A > 30 prior. When we drop the constraint on A^the central 
value of ^ is unexpectedly large, and ^ = is apparently excluded at the 2a level by the 
WMAPII+SDSS dataset. When we compute A^ for these chains we find a distribution 
peaked at 10-15 e- folds, consistent with [l^, but in conflict with successful cosmological 
inflation. We then turn to the pivot dependence of the bounds on the inflationary and 
spectral parameters, and show that in the absence of a prior on A, choosing a different 
pivot scale (comparing constraints at 0.002 Mpc~^ to those at 0.02 Mpc~^, in this case) 



^For instance, if the post-inflationary universe fias an epocli in whicli a{t) ^ t - as miglit liappen if 
tlie end of inflation is followed by the formation of a "frustrated" network of cosmic strings - A'' and the 
inflationary energy scale are largely uncorrelated (see e.g. [^^). This is admittedly an extreme example, 
and we may explore constraints arising from the connection between A'' and the inflationary scale in future 
work. 



- 4 - 



somewhat tightens the bounds on ^ and while shifting Ug significantly to the red. 

In the case with a "sufficient e-folds" prior, the pivot dependence of the constraints is not 
strong. We stress that these results must be interpreted with caution, and make no claims 
to have found strong evidence for a non-zero value of ^, as we have not yet made use of 
the full cosmological dataset. However, these results help us understand the significance of 
the weak preference for a negative running seen in the WMAPII dataset, and give us cause 
for optimism that future data will allow us to tightly constrain the inflationary parameter 
space. 

This paper is arranged as follows. In Section 2 we briefly review the slow roll hierarchy, 
and its use in Monte Carlo Markov Chain analyses. Section 3 contains the constraints on 
the inflationary parameter space derived from our chains. In Section 4 we give our inter- 
pretation and analysis, including a discussion of the pivot-dependence of our constraints 
on the inflationary parameter space. 



2. Imposing A Single-Field Slow- Roll Inflationary Prior 



The dynamics of single field inflation can be written in the Hamilton-Jacobi form, where 
overdots correspond to time derivatives and primes denote derivatives with respect to 4> 
1 19]. The scale dependence of the HSR parameters ^Xh is given by 
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The usual slow roll parameters are rj = ^Xh and ^ = ^Xh- If all terms with £ > M are zero 
at some fiducial point, these differential equations ensure they vanish at all other times. 



Liddle showed that the hierarchy can be solved exactly when truncated at order M |21|, so 
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The Bi are specified by the initial values of the HSR parameters, 
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where the subscript refers to their values at the moment the fiducial mode /cq leaves the 
horizon, and (/> = 0o = 0. We choose = 0.002 Mpc~^. At this point, we emphasize that 
this fiducial scale is not a pivot point - it ties the potential directly to a physical scale in 
the universe, and the closed nature of the truncation of the fiow hierarchy then allows us 
to compute the primordial spectrum at any other physical scale. The potential is 
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while N is given by 

dN _ An H 

and (j) and k are related by 

# _ mpi y/e 
dink 2/^1 -e' ^ ' ' 

Finally, we express the scalar and tensor primordial power spectra in terms of these param- 
eters. A first order expansion around an exact solution for the case of power law inflation 
1 37] gives 



[1 - (2C + l)e + Cr,f ( H ^ 
= 



vre \Tnpi 

2 



(2.8) 

k=aH 



P.= [l-(C + l)ef ^ 



(2.9) 

vr \mpij 

k=aH 

where r/ = U^, C = -2 + ln2 + 7Ps -0.729637 and 7 is the Euler-Mascheroni constant. 



One minor inconsistency in our approach is that (|2.8| ) and (2^) are calculated using the slow 
roll approximation, and truncated at second order. In practice, we do not believe this is a 
significant drawback, as the datasets we are working with provide meaningful constraints 
on (at most) the lowest three slow-roll parameters. In the context of future high-precision 
data, this problem could be surmounted by either using a higher order expansion for the 
spectra, or even integrating the mode equations directly. 

We can write the usual observables in terms of the slow roll parameters, 

n, = 1 + 2r? - 4e - 2(1 + C)e^ - ^(3 - 5C)ery + ^(3 - C)^ (2.10) 
r = 16e [1 2C(e - ??)] (2.11) 

'"•^ ' ^2e + 8.2-10e, + !^eC+^e4 (2-12) 



dink 1- 

nt = -2e - {3 + C)e^ + {1 + C)er] (2.13) 

where C = 4(ln 2-I-7) —5, and we have introduced the customary notation a = dug/dln k; r 
is the tensor /scalar ratio^'^ and nt is the tensor spectral index. We retain all terms in 
up to quadratic order in the slow roll parameters, anticipating that ^ may be as large as e 
or rj. Note that these are derived quantities and are never directly used in the calculation 
of power spectra in this work. 

It is also worth emphasizing that the inflationary single-field consistency condition 
nt = —r/8 only strictly holds at the pivot point when the spectra are specified in terms of 
their respective indices and running. This is easily seen by writing out the scale dependence 
of the standard power-law variables: 

/ I, \ns{ko)-l+^{dns/dlnk)ln{k/ko) 

Pn{k) = Pn{ko) (^^J , (2.14) 

/ h \nt(kQ)+\(dnt/d\nk)\n(k/ko) 

Ph{k) = Phiko) i-^] . (2.15) 



Beware of the distinction between C and C - these quantities are sometimes confused in the literature. 
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Hence, while it is common to assume that drit/dlnk = and nt = —r(ko)/8 is constant, 
the scale-dependence of the tensor/scalar ratio r{k) = Ph{k) / Pfi{k) is given by 

/ t, \ nsiko)-l+^idns/dlnk)ln{k/ko)+r{ko)/8 

r{k) = r(fco) i^-^J (2.16) 

and it is apparent that the consistency condition does not hold away from the pivot. This 
is an unsatisfactory state of affairs if one wishes to impose a single-field slow-roll prior, and 
slow roll reconstruction does not suffer from this shortcoming. 



An alternative to slow roll reconstruction is given by the work of Leach et al. |26, 27 



28, 29 1, who specify the spectrum in terms of the slow roll parameters, but then compute 
its scale dependence via equations 2.10| to 2.12[ Conversely, we dispense with an explicit 
pivot point and compute the spectra directly from equations and |2.9| , while using the 
flow equations to obtain the appropriate scale dependent values of the slow roll parameters. 
When only e and rj are included in the analysis, these two approaches (i.e. ours vs. Leach 
et aZ.'s) overlap to with a few percent over the full observable range. However, when ^ is 
included the discrepancy can be much larger. For instance with a pivot at k = 0.002 Mpc~^, 
the spectra computed using the formalism of [26, 27, 28] differ from that obtained using 
the flow equations by as much as 10% at very short scales, as illustrated in Figure m. This 



discrepancy could be corrected by extending the approach of |g6|, |27|, |2^] to higher order in 
slow roll. 

We perform standard MCMC fits to cosmological data, using modified versions of 
CAMB^i [||] to provide the CMB spectra and COSMOMC^^ [|| to handle the Markov 



^^http : / / camb . info/ 

^^http : / / cosmologist . info/ cosmomc/ 
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Figure 1: The difference between the scalar power spectrum computed using ns(e, 77,^) and the 
full HSR parameters at 0.1 Mpc~^ (left) and 1 Mpc~^ (right), for ^ = 0.01. A positive difference 
corresponds to the HSR result being larger than the ns{e,r],^) value. The "missing" region corre- 
sponds to parameter choices which would lead to less than 15 additional e- folds of inflation. Note 
that ^ = 0.01 is significantly smaller than the central value implied by the data, and the discrepancy 
grows with ^. 
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Chain # 


Dataset Combination 


E-folds Constraint 


1. 


WMAPII +SDSS 




2. 


WMAPII 


N >30 


3. 


WMAPII +SDSS 


N >30 



Table 1: MCMC Parameter fits described in this paper. 



Parameter 


WMAPII+SDSS 


WMAPII 


WMAPII±SDSS 




(no N restriction) 


{N > 30) 


{N > 30) 


e 


< 0.013 (95% CL) 


< 0.033 (95% CL) 


< 0.015 (95% CL) 


V 


0.029 ±0.021 




0.0 ±0.014 




0.029 ±0.013 


0.0 ±0.003 


0.0 ± 0.002 


ln[10i%] 


3.17 ±0.06 


3.03 ±0.10 


3.17 ±0.06 



Table 2: Constraints on HSR parameters defined at k — 0.002 Mpc ^. Constraints are at the 68% 
confidence level unless otherwise noted. 



chains. All our chains include uniform priors on the slow roll parameters {eQ,rjQ, ^q}, where 
eo, rjo, and are the values at ko = 0.002 Mpc^^. The primordial curvature power 
spectrum is normalized at ko by setting 

and the uniform prior is on ln[10^''^s]. Hereafter, we drop the subscripts for brevity. 
The late time parameters are Q^h'^, QcDAih'^, h (the uniform prior is on 9a) and r, where 
these variables have their usual meanings. Since we impose an inflationary prior we take 
^Totai = 1 so the contribution from dark energy is JIa = ^ — ^b—^CDM, and we set w\ = — 1. 
Instantaneous reionization is assumed. To constrain this parameter set, WMAPII and 
SDSS data are used. We use the v2pl version of the WMAP 3 year likelihood code^'^ with 
the same options that were employed in ||2^, and the SDSS [p^ j likelihood routine provided 
in the May 2006 version of COSMOMC. The latter marginalizes over the normalization bias 
of the SDSS galaxy power spectrum with respect to the underlying matter power spectrum. 
We do not impose an age prior, nor do we employ the HST prior for the Hubble constant. 
In each case we ran 8 individual chains. 



3. Results 

We have run three different sets of chains, as summarized in Table |l|. We estimate the 
parameters for WMAPII alone, and for WMAPII±SDSS, with no constraint on N and with 
the requirement that > 30. In each case we tested convergence via the usual Gelman &: 



Rubin 1 40 1 criterion. The chains that drew only on the WMAPII dataset took a very long 
time to converge. Without the > 30 constraint we were not close to convergence even 
after taking over 500,000 accepted steps. For this reason we do not report constraints for 



^http : / /lambda. gsfc .nasa.gov/product/map/dr2/ 
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1-' o v o Tvi ctT ay 
JT di dllic Lei 


(no restriction) 


WM A PTT 

{N > 30) 


vv iVlr\l lln^Ol^OO 

{N > 30) 


Us 


1.08 ±0.05 


0.99 ± 0.02 


0.98 ±0.02 


dns/dln k 




0.004 ± 0.006 


n n+O-004 

'-'•^-0.003 


r 


< 0.22 (95% CL) 


< 0.56 (95% CL) 


< 0.24 (95% CL) 


nt 


> -0.026 (95% CL) 


> -0.067 (95% CL) 


> -0.030 (95% CL) 



Table 3: Constraints on derived primordial parameters defined at fc = 0.002 Mpc ^. Constraints 
are at the 68% confidence level unless otherwise noted. 



the WMAPII dataset on its own, without a prior on N. When SDSS data was included, 
convergence was comparatively prompt, reflecting the greater leverage provided by the 
combination of CMB and LSS data. 

The prior (or lack thereof) on N is implemented as follows. Each draw of {e, rj, ^} in the 
chains maps uniquely onto a potential. For this potential, we can evolve it forward in time 
from when the fiducial scale kQ leaves the horizon and compute when e = 1, corresponding 
to the end of inflation. In the case with no prior on N, we only check that inflation lasts 
for at least the few e-folds encompassing the CMB and LSS data that we actually use - if 
inflation ends within this observable window, that model is discarded. For the case where 
N > 30, we require that the potential so specified provides at least an additional 30 e-folds 
of inflation after the fiducial scale leaves the horizon. Otherwise it is rejected from the 
chain. 

The parameter bounds derived from each set of chains are given in Tables ^ and ^, while 
Figures m through ^ show the results derived from the chains we ran. In particular. Figure ^ 
shows the limits on the slow roll parameters derived from Set 1, where we are fitting to 
the combined power spectrum derived from WMAPII and SDSS, with no constraint on 
the number of e-folds. In Figure |3| we show the same results, after they have been mapped 
into the "standard" slow roll variables. These plots are analogous to those shown in p2| , 
except for the presence of the ^ parameter. Comparing these plots with those from our 
previous paper, we see that including ^ greatly increases the allowed range of values of the 
running, compared to the {e, ry} chains. Indeed, a naive interpretation of the data suggests 
that there is the evidence for ^ > beyond the 2a level - but there are a number of issues 
surrounding this result, and we stress that we make no claim to have found a non-trivial 
value of ^. In Figure Q we compare these results to chains run by Spergel et al. for 
the same dataset. The two sets of chains are not exactly equivalent as the Spergel results 
have an additional constraint on the bias for the SDSS, which we do not include in the 
chains described here. We again see that a naive inspection of the allowed portion of the 
parameter space shows that working with an inflationary prior imposed via the use of the 
slow roll parameters in the chains leads to a stronger preference for a negative running 
than was found by Spergel et al. A further difference is that upper bound on the tensor 
contribution found using a slow roll prior is considerably lower than that obtained from the 
{ra^jT, ^1^^ } parametrization of the spectrum. This phenomenon was also observed in |22| 
in the absence of the running, and appears to be related to the fact that the degeneracy 
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between e and r] in terms of their contribution to Ug is broken by the way these quantities 
contribute to the running index and r. 

A further, more general, message to be drawn from the two sets of chains shown in 
Figure ^ is that the strength of the evidence for a running spectrum is dependent on the way 
we parameterize the running. One might assume that imposing a strong inflationary prior 
by directly including the slow roll parameters in the chains would decrease the available 
region of parameter space, relative to that found from using the conventional {ns,r, -j^^} 
formalism. However, we see that the two parameter volumes are slightly offset, and neither 
one is a strict subset of the other. This reminds us that parameterizing the power spectrum 
in terms of the index and running is simply an empirical fit to what we assume is a weakly 
scale dependent function, but this characterization of the spectrum is not fundamental. 
A further caveat is that we do not impose a bias prior for the SDSS power spectrum in 
our chains, unlike those of Q. This constraint could have a significant impact on the 
estimate of the running, since it effectively fixes the normalization of the SDSS spectrum. 
We cannot distinguish at present between which of these effects - the change in priors and 
parameterization, or the SDSS bias prior - is most responsible for this difference from 




0.01 0.02 0.05 0.1 




0.01 0.02 0.05 0.1 0.020.040.06 




0.01 0.02 0.05 0.1 0.020.040.06 3 3,2 

E Ti ^ log[io™Aj 



Figure 2: Primordial parameter constraints for WMAPII+SDSS (no N constraint). The top plot 
in each column shows the probability distribution function for each of the parameters, while the 
other plots show their joint 68% and 95% confidence levels. 
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This would require re-running our model with the SDSS bias prior; we leave this issue for 
possible future study. 

As noted above, using the slow roll prior puts tighter constraints on the amplitude 
of the tensor spectrum than the {ns,r, chains. If one looks carefully at Figures ^ 

and ^, one sees that e and r are peaked at a non-zero value. This result is not significant, 
but it is worth pointing out that if the potential does have a complicated shape, it may 
be possible to measure e without first detecting the primordial tensor contribution to the 
CMB. Recall that e determines how fast the inflaton field evolves, whereas the scalar index 
is (to lowest order) a linear combination of e and rj. Finally, we remind the reader that 
in the conventional parameterization, the slope and amplitude of the primordial tensor 
spectrum are determined by r. Consequently, the plots in Figure ^ do not contain nj, 
whereas in Figure |2| we can plot a distribution in the {nf,r} plane, since second order 
terms in slow roll ensure that nt is not strictly proportional to r. 

A novel feature of slow roll reconstruction is that we obtain information about the likely 
number of e- folds of inflation, estimated in terms of e, r/ and ^ measured at some fiducial 
moment during inflation. Conversely, we can invert this process and insert a prior on the 
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Figure 3: Primordial parameter constraints for WMAPII+SDSS (no N constraint) in conventional 
power law variables. The top plot in each column shows the probability distribution function for 
each of the parameters, while the other plots show their joint 68% and 95% confidence levels. The 
power law variables are derived parameters which are not used in the MCMC chains themselves. 
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number of e-folds into the parameter estimation algorithm. In practice, we only work with 
lower limits on A^. While it is not possible for an arbitrary amount of inflation to occur 
after cosmological scales leave the horizon, there is a lower bound on the number of e-folds 
if the universe is to successfully reheat. In many models, the end of inflation occurs when 
an instability develops in a direction in field-space that is orthogonal to the inflationary 
trajectory. In that case there is no "advance warning" that inflation is about to end, and 
the slow roll parameters may remain small right up until the transition to non-inflationary 
expansion. Here a naive computation of N using the values the slow roll parameters as 
the CMB scales leave the horizon yields a misleadingly large number. Conversely, a low 
value of N (which is practice means anything less than 30), implies that inflation will end 
too quickly unless some further effect comes into play. This could be a contribution from 
higher order terms in slow roll, or a second period of inflation - both of which are excluded 
by our prior. Consequently, we can usefully include a lower bound on in our fits, but 
we do not impose an upper bound. 

In Sets 2 and 3 we impose the prior > 30 on both WMAPII, and the combination 
WMAPII+SDSS, and plot the results in Figure |5[ These chains put much tighter limits 
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Figure 4: Primordial parameter constraints for WMAPII-I-SDSS (no N constraint) in conventional 
power law variables (black), compared with the most directly comparable results from |^ using 
chains obtained from http://lambda.gsfc.iiasa.gov (blue). The Spergel et al. chains contain an 
additional bias constraint on the SDSS data which our runs did not include. 
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on ^ than Set 1, which has no constraint on A^. Further, the addition of the SDSS data 
to WMAPII significantly tightens the parameter constraints. Comparing our results here 
with Figure 5 of |22|, where we perform a fit to the first two slow roll parameters, we see 
that the permitted region of the {e, rj} plane does not change substantially when the ^ 
parameter is added along with the > 30 constraint. Without the > 30 constraint 
adding ^ significantly expands the allowed range of {e,r]}. 

In Figure ^ the two sets of > 30 chains are converted into constraints on , r and 
dris rpj^g constraints on jrH: are an order of magnitude tighter than those we obtained in 



dink' dink 

the absence of the A^ > 30 prior. The tightest published constraints on fP"", are those given 



by |16], who include Lyman-a and SNla data, along with the information from SDSS and 
WMAPII that we employ here. This combination yields = (— 1.5± 1.2) x 10^^, where 
the error range spans the 68% confidence interval. We present our constraints on the slow 
roll parameters in Table ^ while Table ^ contains the corresponding limits on the spectral 
parameters. With the A^ > 30 constraint, our bound on is 6 x 10~^ - half that of |16|, 
and it shrinks again when we include the SDSS data, and both analyses return a running 
that is peaked very close to zero. Conversely, our constraint on Ug is noticeably broader, 
and shifted toward the blue, relative to that of ||l6| - we shall see later that this may be 
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Figure 5: Primordial parameter constraints for WMAPII+SDSS, (black) and WMAPII (red), 
imposing a prior A'^ > 30. The top plot in each column shows the probability distribution function 
for each of the parameters, while the other plots show their joint 68% and 95% confidence levels. 
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due to the fact that they report their constraints at a pivot point (0.05 Mpc~^) which is 
at significantly smaher scales than ours. We do not include the dark energy equation of 
state (w) in our analysis, whereas this is a free parameter in jl^. The differences between 
our results and those of [^] reflects different underlying assumptions, rather than any 
intrinsic "conflict" between the two analyses. The authors of set out to work with 
as large a range of data sources as possible, whereas our approach is to work with only 
one or two probes of the fundamental power spectrum, and to impose a strong theoretical 
prior - namely slow roll inflation - on the mechanism responsible for the generation of 
perturbations. 



4. Interpretation and Analysis 



While slow roll reconstruction hinges on a very simple idea - inserting the slow roll param- 
eters directly into the cosmological parameter estimation process - interpreting the results 
calls for a good deal of sophistication. Without a constraint on the number of e-folds, we 
see an apparent preference for a large, positive value of ^, which translates into substantial 

-0.05 



and negative values of 



dink' 



Recalling equation 2.14, we see that if 
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Figure 6: Primordial parameter constraints for WMAPIH-SDSS, (black) and WMAPII (red), 
imposing a prior N > 30, in conventional variables. The top plot in each column shows the 
probability distribution function for each of the parameters, while the other plots show their joint 
68% and 95% confidence levels. The power law variables are derived parameters which are not used 
in the MCMC chains themselves. 
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the center of the distribution derived from the combination WMAPII+SDSS), the running 
term dominates this expression as we move away from feo- As is weU known, conventional 
cosmological observations probe the primordial spectrum over a range of around three or- 
ders of magnitude in scale. This may stretch to five orders in the future as we better 
probe the Lyman-a forest or obtain a power spectrum from observations of high redshift 
21cm emission. With |^|^| ~ 0.05, the second term in the exponent of equation 2.14 will 
dominate over the first for k values that fall well inside the cosmological "window", and 
will dramatically modify the power relative to that seen at the pivot, /cq. For example, 
with Us ~ 1, ln(fe//co) ~ 10 and —0.05, the power spectrum at fc = e^'^ko is an 

order of magnitude smaller than at ko. However, if we have a spectrum that is strongly 
deficient in short scale power, the other variables in our parameter set will move to offset 
this deficit. As we can see from Figure this is exactly what happens here, as has 
increased relative to its concordance value, where fl\ and Hq have moved downwards. The 
former is simply a result of our assumption of a fiat universe, so that Jlxotai = 1, and the 
increased Om compensates for the reduced power in the perturbation spectrum at short 
scales. 

Imposing additional constraints on or Q\ would put pressure on the models that 
produce large and small Q\. In Figure § we show the distribution of {n^, dn^/dln fc} 
for Set 1 (WMAPII+SDSS), with the points color-coded by Qru- This plot illustrates the 
correlation between an extreme running and an anomalously large Consequently, a 
strong upper bound on Jim, from observations independent of the CMB and LSS data we 
employed in our MCMC analysis would exclude the models with the most extreme running. 
Consequently, the long tail of the allowed region in the {n^, dus/dlnk} plane for Set 1 can 
be understood as a marginalization artefact. 

In models with large ^, rj and e are typically strongly scale-dependent, thanks to 
feedback from the ^ term in the flow equations. Once e is equal to unity, inflation will 
come to an end. As we have explained, imposing > 30 leads to a dramatic change 
in the permitted region of parameter space, consistent with our qualitative analysis in 
1 18]. In Figure ^, we show the remaining number of e-folds as a function ^, obtained 
by post-processing our chains. Note that only models in which inflation explicitly ends 
are retained from the chains in order to make these plots. Models which are drawn to 
a "late-time attractor" and inflate forever (thus needing an orthogonal mechanism to end 
inflation) cannot be assigned an unambiguous number of remaining e-folds after the flducial 
scale leaves the horizon. For this reason they are excluded from the figure (and we have 
checked that the Markov chains cut in this way still retain their convergence properties). 
We immediately see that Set 1 - chains which ran without any constraint on - has a 
distribution in the {^, A^} plane that peaks near A^ = 10, and excludes A^ = 30 at more 
than 2a. On the other hand, the corresponding distribution for chains with A^ > 30 is 
peaked around A^ = 30, but A^ = 60 is no longer excluded. 

While we have dubbed the methodology employed here "slow roll reconstruction" , we 
are not in a position to present an unambiguous form for the infiationary potential. With- 
out a non-zero lower bound on e the height of the potential is essentially undetermined. 
Secondly, as we noted above, our chains prefer a positive value of ^, but we emphasize 
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that this is stih weh short of a convincing detection of a running, nor have we performed 
a Bayesian analysis (e.g., [41, |4^) to test the significance of the inclusion of the ^ param- 
eter. On the other hand, the underlying power spectrum is uniquely determined by any 



given choice of {As,e,rj,(^}, and in Figure 10 we plot the permitted range of P{k) for the 



three sets of chains. In particular, without the > 30 prior we see the preference for 
a scale dependent spectral index in the WMAP+SDSS dataset. Conversely, it is easy to 
understand why including the Lyman-a dataset has such a dramatic impact on the allow- 
able running, since this samples the power spectrum at very short scales, and thus has 
substantial "leverage" . 

Figures |9| and |l^ show that large values of ^ have a strong tendency to make inflation 
end after an insufficient number of e-folds to satisfy cosmological requirements. These 
figures imply that in order to satisfy the requirement that > 30, the allowed range 
of ^ from the data becomes dramatically restricted, leading to a small allowed range on 



dink 

While we specify the initial values of the slow roll parameters at a particular scale, 
we are able to post-process our chains to produce plots for the slow roll and spectral 
parameters at any arbitrary pivot (of course ensuring that the convergence criteria for the 
post-processed chains remain satisfied). This is of particular importance when deal with 
a running index, since is an explicit function of k, and the bounds on thus depend 
explicitly on k^. In our chains we took ko = 0.002 Mpc"\ but we now consider how the 
bounds on the parameter space vary as we translate the pivot. 



It is apparent from Figure 11 that modifying kg significantly changes the constraints 

dns 
dink' 



on the slow roll parameters and and . Here we have chosen a pivot near the "waist" 



of the P{k) distribution in Figure or fco = 0.02 Mpc"^ For the case where no 
prior was imposed, we now find that there is a considerable shift to the red in Ug but the 
bounds on ^ and -^f^ are not that much tighter than with ko = 0.002 Mpc^^. However, 
the Ug — degeneracy has rotated somewhat, encapsulating both the tendency seen in 



the top panel of Figure for the power spectrum to run from blue on large scales to 
red on small scales, and the fact that specifying the constraints at the "waist" decorrelate 
these two parameters as much as possible. When we impose the prior, the constraints 
do not change significantly when we translate the pivot - in this case the running was 
already tightly constrained and the pivot dependence of is correspondingly small. We 
can compare these plots directly to specific inflationary models, and we have superimposed 



the trajectories for several canonical models on the plots in Figure 11. A representative 



natural inflation [43| model and m?(jp' inflation lie inside the 95% confidence level with 
and without the A^ prior for the combination of WMAPII-I-SDSS. Conversely, Xcf/^ lies 
outside this confidence level for both cases. We do not show results for power-law or 
hybrid inflation. In the former case, Hg and the slow roll parameters do not evolve with 
time, and so any given model would be represented by a single point on these plots, rather 
than a trajectory. Likewise, in hybrid models the end of inflation is not specified by the 



Natural inflation contains a free parameter / that appears in the denominator of the cosine potential. 
Bounds derived from the WMAPI dataset on natural inflation are discussed in 144|. 



-16- 



value of the slow roll parameters, so we cannot unambiguously plot their trajectories as a 
function of N. 

Looking at Figure ^ we are again reminded that interpreting MCMC estimates of 
cosmological parameters is a complicated task. At present, the publicly available Lyman- 
a likelihood codes are not compatible with the Hubble slow roll formalism used here. 
Since these data probe the fundamental spectrum at small scales, incorporating it into 
our analysis would greatly enhance our ability to constrain the spectrum - particularly for 
models which are strongly scale dependent. Likewise, we expect that the release of high 
quality CMB spectra for i > 1000 will greatly enhance the ability of slow roll reconstruction 
to constrain the inflationary parameter space. 



A further lesson to be drawn from Figure 11 is that the data is now approaching 
the point where we must directly address the dependence of the parameters in simple 
inflationary models. In |4^], Kinney and Riotto show that our ignorance of reheating 
physics means that we cannot make a unique mapping between the remaining number of 
e- folds A^ and a given comoving scale k. Looking at the trajectories we plot, we see that the 
N dependence of the slow roll and spectral parameters over a ten e-fold range is only a few 
times smaller than the widths of the parameter distributions we recover from our MCMC 
analyses. The quality of cosmological data is thus reaching the point where this degeneracy 
will need to be included explicitly in any experimental tests of specific inflationary models 
~ especially those where the field point evolves comparatively rapidly. 

5. Discussion 

In this paper we have used slow roll reconstruction to extract bounds on the inflationary 
potential from observational data. These bounds make strong use of an inflationary prior, 
and thus apply to cosmological models where the density perturbation spectrum is laid 
down by single, slow rolling, minimally coupled inflaton field. Our results align with the 
broad picture supplied by previous analyses of the WMAPII dataset: there is noticeable 
but by no means compelling evidence for a significant scale dependence (running) in the 
power spectrum when one considers WMAPII and SDSS on their own. However, adding 
a prior constraint that rules out models with an unacceptably small number of e-folds 
significantly tightens the bounds on any possible running. 

The combination of WMAPII and SDSS on their own are compatible with a broad 
class of perturbation spectra, but this range shrinks substantially when further constraints 
are added to the estimation process. In particular, we consider inflationary models where 
at least 30 e-folds of inflation occur after CMB scales have left their horizon. Similar results 
are seen in conventional analyses of the spectrum when Lyman-a forest data is included 



16]. We further analyze the impact on the choice of pivot on the quoted ranges for the 
spectral parameters. Slow roll reconstruction provides an estimate for the inflationary 
potential over its entire range, thanks to the closed nature of the truncated slow roll 
hierarchy. Consequently, we actually have access to the underlying power spectrum and 
inflationary potential, although these quantities are fully specified by the truncated slow roll 
hierarchy - which here consists of three independent parameters. However, we can use the 
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flow equations to post-process our chains to give the constraints on the standard spectral 
parameters at an arbitrary pivot. We make use of this feature, and study the change in 
constraints on the inflationary parameter space when changing kg from 0.002 Mpc~^ to 
0.02 Mpc"^ Without a prior on the number of e-folds A^, the constraints on Ug and 
are found to be significantly dependent on the pivot scale. The large running seen without 
this prior cannot be removed by shifting the pivot - changing ko shifts Ug well to the red, 
so the hints of an overall scale dependence cannot be expunged simply by changing ko. 

One lesson from this analysis is that pivot-dependent estimates of power-law observ- 
ables, while valid, need to be interpreted with great care. For example, in the case of the 
constraints presented at A: = 0.02 Mpc"^ for WMAPII-I-SDSS without any prior on N, the 
fact that this pivot scale is an inflection point "hides" a tendency to run from blue on large 
scales to red on small scales as you move away from this pivot. That is, when there is a 
strong running, any slice in k of the primordial power spectra can look quite different in n^, 
r, space. In our view, it is a better aid to interpretation to present constraints over 



the entire k range implied by one's parameterization, as in Figure 10 ~ in principle, this is 



a simple thing to do in the MCMC methodology for any parameterization of P{k) - and 
then compare with specific infiationary models over this entire range of scales. What does 
matter is whether the constraints at a wide range of k either encompass or exclude the pre- 
diction from that specific potential. For absolute completeness, the "band" of constraints 
from data on P{k) should be compared with a "band" for the model under consideration, 
encompassing the theoretical uncertainty in both the model parameters and the reheating 
process [|45|| . 

Where they are directly comparable, our results broadly agree with those of Finelli et 
al. [^9|. However, the authors of use a different formulation of the slow roll expansion 
from the one adopted here, which allows their ea to be of order unity, and they find 
they must pay careful attention to their prior for this parameter. Conversely, our "third" 
parameter, ^ needs no such special handling, and we simply take all our slow roll parameters 
to lie within [—1, 1], which is far larger than the range allowed by the data. Secondly, our 
running (a) is effectively scale-dependent, since ^ is a function of scale, but [29| specify 
their running at a fixed pivot, and it is not a function of scale. Finally, by using slow roll 
reconstruction we are able to include constraints on the remaining number of e-folds A^, a 
possibility which is not considered in the analysis of |^^. The tension between the running 
and the number of e-folds confirms our theoretical analysis in [^] . 

Algebraically, the running index is dominated by ^, since e and r] are both required 
by the data to be relatively small, and these terms only contribute quadratically to ^j^. 
When this parameter becomes large and positive the fiow equations ensure that e and r/ 
will grow quickly. In this case e becomes equal to unity in ~ 15 e-folds, and infiation 
terminates far too quickly for it to solve the usual cosmological problems. In this analysis 
we have post-processed our chains to obtain likelihood contours for N in the (^, N) plane, 
and we see that does in fact peak around 10 e-folds. Conversely, if we impose a prior on 
N ^ the permitted range of -^-^ shrinks considerably. 

^^In the formalism of £3 only appears with a small factor, so while this parameter need not be small, 
its overall contribution is still typically sub-dominant. 
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At large values, ^ has a long degeneracy because of the lack of constraining power in 
the current dataset. This is possibly a marginalization artifact, as we found that a large 
running is correlated with values of Qm beyond the range seen in the standard concordance 
cosmology. If "simple" models of inflation are a good description of the universe, the rising 
quality of astrophysical data will eventually break this degeneracy, and the running will fall 
somewhere within the range found when the > 30 prior is applied here. Conversely if 
the data contracts around the large negative median value of the running, we would learn 
that inflation is non- minimal in some way, as discussed in [18|. In either case, slow roll 
reconstruction will be able to put tighter constraints on N as the constraints on ^ tighten. 

In a subsequent paper we plan to assess the ability slow roll reconstruction to put 
bounds on the inflationary parameter space with the data from different proposed and pro- 
jected experiments. In addition, we will need to include A3, the fourth slow roll parameter, 
in these calculations, to be sure than any conclusions we reach are not a function of a 
premature truncation of the slow roll hierarchy. On the other hand, if this latter analysis 
shows that A3 is necessary to describe the data four independent parameters would be 
needed to describe the inflationary potential, which may be an indication that inflation is 
non-minimal. 

Finally we stress that while we find evidence that a significant breaking of scale- 
invariance in the primordial power spectrum is consistent with the 3-year WMAP dataset 
in the presence of an inflationary prior, there is no compelling evidence for a running 
spectral index. However, the analysis here gives us cause for optimism that this question 
can be settled in the near future, and that slow roll reconstruction provides an elegant and 
powerful framework for analyzing the cosmological constraints on slow roll inflation. 
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Figure 7: Marginalized ID posterior distributions for the "late time variables" for 
WMAPII+SDSS, imposing a prior TV > 30 (dashed) and with no constraint on N (solid), r 
and as are remarkably stable to the addition of the prior; most other parameters experience shifts 
in the mean at the level of a u or more. The greatest effect is on Sl^ (see the discussion in text for 
details). 
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Figure 8: The joint 68% and 95% confidence levels in the Us vs dus/dlnk plane derived from the 
WMAPII+SDSS chains with no N constraint, showing fim = l^b + ^Icdm, the late time variable 
which is most correlated with this eigenvector. This gives an indication of why the SDSS data does 
not help to break this degeneracy more strongly, despite the greater leverage on fc-coverage one gets 
on adding it. 
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Figure 9: The joint 68% and 95% confidence levels on the number of e-folds after the scale 
ko = 0.002 Mpc"^ leaves the horizon tiU the end of inflation vs ^ for WMAPII+SDSS (black) and 
WMAPII (red); no constraint on N (left), > 30 (right). Only models from each set of chains 
where inflation explicitly ends are included. 
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Figure 10: Reconstructed 68% (shaded dark) and 95% (shaded light) constraints on the power 
spectrum of the primordial curvature perturbation as a function of fc, shown for WMAPII+SDSS 
with no N constraint (top), WMAPII with an > 30 prior (bottom left), and WMAPII+SDSS 
with an A^ > 30 prior (bottom right). The solid line denotes the mean. The scales spanned by 
WMAP, SDSS and Lyman a data (not used in this analysis) are shown. The results with no A^ 
constraint show a tendency of the power spectrum to run from blue at large scales to red at small 
scales, which disappears when the prior on A^ is added. The dashed line shows the 95% CL upper 
limit on the power spectrum of the primordial tensor perturbation. 
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Figure 11: We present bounds on the slow roll variables and spectral parameters, with the pivot 
moved to ko = 0.02 Mpc~^. The constraints are drawn from our WMAPII+SDSS chains, The blue 
regions are derived from Set 1 (no prior on N), and the brown from Set 3 {N > 30). We show 
68% and 95% confidence intervals. We superimpose the "trajectories" for three generic slow roll 
models, (dotted), m'^cj)'^ (soHd), and a representative natural inflation model (dashed), showing 
how the predictions of these models depend on iV: N = 30 (triangle), = 40 (square), iV = 50 
(star), A'^ = 60 (circle). The figures containing running and ^ only show the m^0^ model, as the 
other models make similar predictions for the running. 
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